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^^ ■ Abstract 

Let U be a multiply-connected region in R 2 with smooth boundary. Let P £ be a polyomino in eZ 2 
approximating U as e — » 0. We show that, for certain boundary conditions on P,,,the height distribution 
on a random domino tiling (dimer covering) of P e is conformally invariant in the limit as e tends to 0, 
in the sense that the distribution of heights of boundary components (or rather, the difference of the 
heights from their mean values) only depends on the conformal type of U. The mean height is not strictly 
conformally invariant but transforms analytically under conformal mappings in a simple way. The mean 
height and all the moments are explicitly evaluated. 
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1 Introduction 

Conformal invariance of a lattice-based statistical mechanical system is a symmetry property of the 
system at large scales. It says that, in the limit as the lattice spacing e tends to 0, macroscopic quantities 
associated with the system transform covariantly under conformal maps of the domain. 

Conformal invariance for statistical mechanical lattice models is a physical principle which until now 
has not been proved except in certain models which were tailored to be conformally invariant M (recently 
in |2j Benjamini and Schramm prove conformal invariance in a discrete, but non-lattice, percolation 
model). Nonetheless conformal invariance is an extremely powerful principle: in the plane, conformally 
invariant models are classified, in a sense, by representations of the Virasoro algebra |Dj . Physicists have 
used this theory fruitfully to compute exact "critical exponents" and other physical quantities associated 
to critical lattice models M . For example, the cycle in Figure [U is believed to have Hausdorff dimension | 
in the limit (see e.g. (l5|) and the path in Figure w is believed to have dimension | Jill]. Although many 
well-known models are believed to be conformally invariant at their critical point, no rigorous techniques 
were known to prove conformal invariance in these models. 

In this paper we deal with the two-dimensional lattice dimer model, or domino tiling model (a domino 
tiling is a tiling with 2x1 and 1x2 rectangles) . We prove that in the limit as the lattice spacing e tends 
to zero, certain macroscopic properties of the tiling are conformally invariant. 

The height function h o n a domino tilin g is an integer-valued function on the vertices in a tiling. 
It is defined below in section |2.jj : see also M, [19| . One can think of a domino tiling of U as a map h 
from U to Z, where for each unit lattice square, the images of the four vertices under h are 4 consecutive 
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integers v, v + 1, v + 2, v + 3. Furthermore each boundary edge of U must have image of length 1 and not 
3. The map h defines and is defined by the tiling: the edges crossed by a domino are those whose image 
under h has length 3. Our main result is the conformal invariance of h for a random tiling: 

Theorem 1 Let U be a bounded, multiply connected domain in C — R 2 with k + 1 smooth boundary 
components, each with a marked point do,di,... ,dk- Let {P E }e>o be a sequence of polyominos, with 
P c G el? , approximating U as described in section 5.S. Let of be a vertex of P e within O(e) of dj. Let 



H c be the uniform measure on domino tilings of P e . Then the joint distribution of the height variations of 
the points er- (that is, the difference of the heights from their mean value) tends to a finite limit which 
is conformally invariant. 

By conformal invariance we mean, if / : U — * U' is a conformal isomorphism then the distribution of 
the height variations of /(dj) is the same as the distribution of the height variations of the dj themselves. 

The mean height of a point of P c is not strictly conformally invariant in the limit: there is an extra 
term coming from the heights on the boundary (Theorem K3). We prove there that the limiting mean 
height is a harmonic function on U whose boundary values depend on the tangent direction of the 
boundary. 

The picture of the height function is completed by understanding the distribution of heights at interior 
points of U. For an interior point x of P e , Theorem |2| below and J13J show that the height h(x) tends 
to a Gaussian with variance clog(i) for a constant c (which can be shown to be \ by a computation 
similar to that in |13|). See below. This variance diverges as e — » 0. On the other hand the proof of 
Theorem hi shows that the moments 



E((/i(xi) - h(xi)){h{x2) - h{x 2 )) ■ ■ ■ (h(x m ) - h(x m ))) 

for distinct Xi tend to a finite and conformally invariant limit. 

Theorem hi can be extended to regions U with piecewise smooth boundary, on condition that at each 
corner the boundary tangents have one-sided limits. See below. 

Figure hi illustrates one consequence of Theorem hi In that figure we took two random domino tilings 
of an annular region (a square with a square hole) . A domino tiling corresponds to a dimer covering, or 
perfect matching, of the underlying graph (a perfect matching is a collection of edges covering each vertex 
exactly once). Two perfect matchings form a union of closed cycles and doubled edges in the graph. One 
can ask about the distribution of the number of cycles separating the inner and outer boundaries of the 
annulus (there is just one such cycle in the figure). The argument of |13| shows that the distribution 
of the height difference between two boundary components for a single domino tiling is directly related 
to the distribution of the number of cycles separating those two components in a union of two tilings. 
Indeed, the expected number of cycles is A times the variance of the height difference. Theorem hi 
therefore implies that the distribution of the number of cycles separating the boundary components from 
each other is conformally invariant. 

Another interpretation of the height function uses the connection between domino tilings and spanning 
trees on Z 2 H. In section hj we relate the height function to the "winding number" of arcs in the 
corresponding spanning tree. 

Theorem |l| follows from a more fundamental result. The coupling function on P E is a function 
C : P e x P f — * C which determines the measure /i e (the uniform measure on the set of all tilings of P € ) in 
the sense that subdeterminants of the coupling function matrix give probabilities of finite configurations 
of dominos occurring in a tiling [h3J. The coupling function is closely related to the Green's function. 
The following is a loose statement of the result. 
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Figure 1: A cycle in a union of two random domino tilings of an annulus. 
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Theorem 2 Let U and {P e } e >o be defined as in TheoremYh. Let v 7^ w be points in the interior of U 
and v^ e ',w^' vertices of P t within 0{e) ofv,w respectively. The coupling function C for domino tilings 
of P e satisfies 

C(v (e) ,W (e) ) = eF j (v,w) + o{e), 

where j = or 1 depending on a parity condition, where Fo and F\ are analytic in the second variable 
and depend only on the conformal type of U . 

For a precise statement see Theorem |l3|. This result has an immediate corollary regarding densities 
of local configurations. 

Corollary 3 In a random tiling of P € , the expected density of occurrences of a local configuration E of 
dominos at a point v in the interior of U is of the form c(E) + cWe(v) + o(e), where c(E) equals the 
density of E in a random tiling of the whole plane el? , and We is a function depending only on the 
conformal type of U . 

The proofs of the above results are given for polyominos with somewhat special boundary conditions. 
We discuss in section H alternate boundary conditions for which it may be possible, using similar methods, 
to prove similar results. We remark that certain restrictions on the boundary are definitely necessary, 
however: in M Cohn, Kenyon and Propp compute the mean height when the height function on the 
boundary is of order ~ . In this case the mean height satisfies a much more complicated non-linear elliptic 
PDE and does not appear to have any simple conformal invariance properties. 

The paper is organized as follows. In section H we define the polyominos, graphs and notations we 
will be using. We also define the height function. In section ra we define discrete analytic functions, and 
show that the coupling function is one. In section H we discuss boundary values of the coupling function. 
In section M we prove Theorem ti. In section H we prove Theorem hi using Theorem H. In section r>/A we 
compute explicitly the average height function on a region. In section m we discuss the connection with 
spanning trees, and in section H we discuss other boundary conditions and give some concluding remarks. 

Acknowledgements. I would like to thank Oded Schramm for many helpful ideas, and the referee for 
several simplifications in section 0. 

2 Definitions 

2.1 Polyominos and their dual graphs 

Let T be the checkerboard tiling of R 2 with unit squares, each square centered at a lattice point of Z 2 , 
and where the square centered at the origin is white. Let Wo be the set of white squares both of whose 
coordinates (the coordinates of the center of the square) are even; let W\ be the set of white squares 
both of whose coordinates are odd. Let Bo be the set of black squares whose coordinates are (1,0) mod 2 
and B\ the set of black squares whose coordinates are (0, 1) mod 2. 

A polyomino is a firiitojj union of unit squares of T bounded by disjoint simple closed lattice paths. 
A corner of (the boundary of) a polyomino is convex if the interior angle is 7r/2; a corner is concave 
if the interior angle is 37r/2. In either case the corner lattice square is the lattice square adjacent to 
the corner, which contains the angle bisector of interior angle. An even polyomino is a polyomino P 
in which all corner squares are of type Bi . Note that this implies that any boundary edge of P whose 
two corners are both convex or both concave has odd length; any boundary edge of P with a convex 
and a concave corner has even length. A polyomino is simply-connected if it has only one boundary 
component. 



Later we will consider some special infinite polyominos. 
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Lemma 4 A simply- connected even polyomino contains one more black square than white square. 

Proof. This is easily proved by induction on the number of corners, starting from the case of a 
rectangle. □ 

A Temperleyan polyomino is a polyomino which is obtained from an even polyomino P as follows. 
Remove from P a black lattice square do adjacent to an edge or corner of the outer boundary of P. For 
each interior boundary component Dj of P, add a black lattice square dj adjacent to an edge of that 
boundary. We assume that dj only borders on a single square of P. See Figure 0. These added squares 




Figure 2: A Temperleyan polyomino. The black squares are in B\, the gray are in Bq. 

will be called exposed squares. Note that do must be in B\ and dj must be in Bo for j > 0. From 
the lemma it follows that a Temperleyan polyomino, even if not simply connected, contains the same 
number of black squares as white squares. 

Let P be an even polyomino, and let Bi(P) be the graph whose vertices are the squares B\ in P, with 
edges connecting all squares at distance 2. Then to each horizontal edge of Bi(P) corresponds a square 
W\ of P (the square it crosses) and to each vertical edge of Bi(P) corresponds a square of type Wo of 
P. To each face of Bi(P) which is not a boundary component of P corresponds a square of P of type 
Bo- The planar graph Bi(P) has a planar dual Bo(P), whose vertices are faces of Bo(P) (squares of type 
Bo), as well as a vertex for each boundary component of P. For a Temperleyan polyomino constructed 
from P, we can still associate the same graphs Bi(P) and Bo(P), but we mark the special vertex do of 
Bi(P) and mark in Bo(P) the special edges adjacent to the di for i > 0. 

Temperley |l7J] gave a bijection between spanning trees on an m x n grid and domino tilings of a 
(2m — 1) x (2n — 1) polyomino with a corner removed. A Temperleyan polyomino is a polyomino which 
arises from a subgraph of the grid by a generalization of his construction, as above, where Bi(P) is the 
subgraph one starts with (see [M). 

The interior dual graph M of a Temperleyan polyomino P is the graph with a vertex for each 
lattice square in P, with edges joining pairs of vertices whose corresponding squares are at distance 1 (in 
other words, it is the dual graph without the boundary vertices). Domino tilings of P are in bijection 
with perfect matchings of its interior dual graph (a perfect matching of a graph is a set of edges such that 
each vertex is an endpoint of exactly one edge). The exposed squares of P are called exposed vertices 
of M. 

The interior dual graph M of a polyomino P is a subgraph of Z 2 and its vertices inherit a coloring 
from the checkerboard coloring of the lattice squares: (x, y) is in Wo if and only if (x, y) = (0, 0) mod 2 
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and so on. We will usually denote a vertex (x, y) £ Z 2 by the complex number x + iy. 

2.2 The height function 

Thurston |19] defines the height function on a domino tiling as follows. The height function is a Z- valued 
function on the vertices of the tiling, defined only up to an additive constant. Start at an arbitrary 
vertex of some domino and define the height there to be 0. For every other vertex v in the tiling, take 
an edge-path 7 from vq to v which follows the boundaries of the dominos. The height along 7 changes 
by ±1 along each edge of 7: if the edge traversed has a black square on its left (which may be exterior 
to the region) then the height increases by 1; if it has a white square on its left then it decreases by 1. 
This defines a height at v. If the tiled region is simply connected, the height is independent of the choice 
of 7 since the height change going around a domino is 0. If the tiled region is not simply connected the 
height is still well-defined as long as each hole contains the same number of black and white squares [J19J] . 
See Figure || 
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Figure 3: Heights in a domino tiling. 

Let M be the interior dual graph of a Temperleyan polyomino P, and take a perfect matching of M. 
A height function on the tiling determines a height function defined on the (non-boundary) faces of M . 
The height function may be defined by assigning an arbitrary value to some face and then applying the 
following rules: for each unmatched edge of M, when following the edge from its black vertex to its white 
vertex, the height of the face on the left minus the height of the face on its right is 1. For matched edges 
this difference is —3. 
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2.2.1 Heights of boundary components 

Let P be a Temperleyan polyomino with boundary components Do , ■ • ■ , Dk where Do is the outer com- 
ponent. Since each Dj encloses the same number of black squares as white squares the net height change 
around each Dj is zero, so the height is well-defined for any tiling of P. 

Given a tiling of P the height function along Dj depends only on the height of any single point on 
Dj. That is, given two points xo, x\ of Dj, let 7 be the path running along Dj from xq to Xi. The height 
difference h(xi) — h(xo) is independent of the tiling since 7 crosses no dominos. Since the height of Dj 
depends only on a single integer value, it makes sense to talk about the height of Dj as a single Z- valued 
random variable. 
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Note how the height changes as you go around a boundary component with the interior of P on your 
left (see Figure ph. Along a straight edge the height alternates between two successive values. Except at 
the exposed vertex, after a right turn the alternating pair decreases by 1, and after a left turn it increases 
by 1 (this follows since all corners are black). This means that for two points on the same boundary 
component, their height is related in a simple way to the amount of winding of the boundary component 
between them (i.e. the number of left turns minus the number of right turns). 

2.3 Tilability of big Temper leyan polyominos 

The Temperleyan polyominos we will be using are those with small lattice spacing which approximate a 
region U with smooth boundary (or piecewise smooth with one-sided limits of tangents at each corner). 
Tilability of such a polyomino can be shown using the following result of Fournier. 



Proposition 5 (|10|) A simply-connected polyomino with the same number of black and white squares 
can be domino-tiled unless there are two boundary vertices x,y whose distance in the ^-metric (length 
of the shortest lattice path from x to y in P) is less than their height difference. 

Actually Fournier's condition is stronger than this (he uses a modified metric) but this will suffice for our 
needs. Also, Fournier only considered simply-connected regions but his argument generalizes to regions 
with many boundary components, as long as a height has been assigned to each component (and one is 
interested in tilings whose height function extends the function already defined on the boundary) . 

Since the region U has a piecewise smooth boundary as defined above, the winding number of the 
boundary path between two points on the same boundary component of U is bounded. As a consequence 
if P e is a Temperleyan polyomino in eZ 2 approximating U (and if locally the boundary of P e follows 
that of U in the sense that they are always directed into the same approximate quadrant), the height 
difference between two points on the same boundary component of P e is approximately the same as the 
winding number of the boundary of U between those two points. Therefore the height function on the 
boundary of P e varies by at most a constant. 

In particular if e is sufficiently small Proposition H and Lemma H show that P e is tilable. 

A more elementary proof of tilability using spanning trees is sketched in section bj 

3 Discrete analytic functions 

The important discrete functions appearing in this article are examples of discrete analytic functions 
(also called monodiffric functions), see fl9[. This section reviews the relevant definitions. Our definition 
is slightly different from the classical definition in M but is equivalent. 

3.1 The dj operator 

We define several operators on Z 2 . The operator d x : C z — » C z is defined by: 

«./(«) = /(« + 1) -/(«- 1). 

Similarly define 

d y f(v) = f{v + i)-f(v-i). 

We define operators 

d z — d x — id y , 

and 

b\ = d x +id y . 
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These operators restrict to operators from C B to C : if / G C B , that is, if / is zero on white vertices, 
then dxf,d y f G C w . Similarly d x ,d y map C w to C s . A discrete analytic function is a function 
F G C B which is real on Bo and pure imaginary on Bi and satisfies b\F — 0. If F — f + ig where 
/ G R s ° and g G R Sl , then F being discrete analytic is equivalent to / and g satisfying the discrete 
Cauchy-Riemann equations 

d x f(v) = d y g(v) for v£ W (1) 

dyf(v) = -d x g(v) for v e Wi. (2) 



(Note that when / G R s ° and g G R Sl , we have d x f, d y g G R w ° and d y f, d x g G 

The function / is called the real part of / + ig, and g is called the imaginary part of / + ig. 

If / + ig satisfies the discrete CR-equations at all but a finite number of (white) vertices, we say that 

f + ig is discrete analytic with poles at those vertices. 

The operators d x ,d v ,d z ,& z - restrict to operators on subgraphs M of Z 2 in a natural way: we consider 

C M to be the subset of C z which consists of functions zero outside of M. We apply the operator and 

then project back to C M . 

3.2 Laplacian 

A simple calculation shows that, if / G R s °, then d z o\f G R s ° and —d z d z f is the Laplacian of / on the 
graph B (Z 2 ). That is, 

-d,dwf{v) = Af(v) = 4f(v) - f(v + 2) - f(v + 2t) - f(v - 2) - f(v - 2i). 

Note that this is 4 times the usual Laplacian since we left out factors of \ in the definition of b\ and 
d z . Often when discussing the discrete Laplacian there is a disagreement about the choice of sign. Here 
we chose the positive (semi-)definite Laplacian, which corresponds in the continuous limit to ~-§-^i — ■§-?■ 

In a similar fashion if g G R Sl then —d z b\g is the Laplacian of g on the graph Bi(Z 2 ). 

In particular if / + ig is discrete analytic on 1? we have d z b\{f + ig) — d z (0) = and so A/ = and 
Ag = 0, where the first A is the Laplacian on Bo(Z 2 ) and the second is the Laplacian on Bi(Z 2 ). 

For a discussion of the boundary behavior of the Laplacian on Bo(-P), see section |Ll| . 

3.3 Weighting the graph 

An alternative way to define discrete analytic functions, which relates more closely with domino tilings, 
is as follows. On the graph 1? put weights on the edges: at each white vertex the four edge weights 
going counterclockwise from the right-going edge are 1, i, —1, —i respectively. See Fig. W. 

Now for a pair of real- valued functions / G R s ° and g G R Sl , the function / + ig is discrete analytic 
if and only if it satisfies K (/ + ig) — 0, where K is the adjacency matrix of 7? with these weights. The 
matrix K is called the Kasteleyn matrix of 1? . Kasteleyn proved that for a finite region the absolute 
value of the determinant of the Kasteleyn matrix is the square of the number of perfect matchings. 
(Usually the Kasteleyn matrix is defined with different weights |12j| : but in fact any choice of complex 
weights of modulus 1 satisfying ac — —bd for the four weights a, b, c, d around a square gives rise to a 
Kasteleyn-like matrix whose determinant counts tilings.) 

When considered as an operator on C B , the operator K is the operator b\. When considered as an 
operator on C w , however, it is — b\. Let K* be the Hermitian conjugate of K. Then the operator K* K 
is acting as the Laplacian on both Bo and Bi. 

Lemma 6 A discrete analytic function on a simply connected Temperleyan region P is determined up 
to an additive (imaginary) constant by its real part. 
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Figure 4: Weights of the Kasteleyn matrix. 



Proof. Note first that Bi(P) is connected. Let / £ R s ° be harmonic on Bo(P). Given the value 
of the imaginary part g at one vertex v £ Pi, the value g(w) for any other vertex w in Pi is uniquely 
determined as follows. Take a path in Bi(P) from v to w. Each edge of the path crosses an edge of 
Bo(P). One of the Cauchy-Riemann equations ((hi) or (H)) at the crossing point determines the difference 
in values of g at the endpoints of this edge. The value g(w) is obtained by summing this difference along 
the path. The harmonicity of / implies that the value g(w) obtained is independent of the path chosen. 
□ 



When the region is not simply connected, in general the conjugate function of a harmonic function 
/ £ R s ° is not single-valued: the "integral" in the above lemma along a path surrounding a hole may 
not be zero. 
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Let M be the interior dual graph of a Temperleyan polyomino P. Let K be the corresponding Kasteleyn 
matrix and let E be a finite collection of disjoint edges of M. Let foi , . . . , b^ and wi, . . . ,Wk be the black 
vertices (respectively white vertices) covered by E. Let jj, be the uniform probability measure on perfect 
matchings of M. 

Theorem 7 (||13]) The ^-probability that E occurs in a perfect matching is given by \ detlK^ 1 )], where 
K^, is the submatrix of K~ x whose rows are indexed bybi,... ,bh and columns are indexed bywi,... ,Wk- 
More precisely, the probability is (— 1)^ Pi+qi aE det(K^ 1 )c, where pt,qi is the index ofbi, resp. Wi, in a 
fixed ordering of the vertices, c — ±1 is a constant depending only on that ordering, and aE is the product 
of the edge weights of the edges E. 

Thus the /^-measures of cylinder sets for perfect matchings on M are determined by this function 
K~ x : M X M — ¥ C, called the coupling function. For historical reasons we denote the coupling 
function with a C. 

Actually this theorem holds for arbitrary bipartite planar graphs, not just those arising from the 
square grid: see M. 

In all of our applications of this theorem we will use only a small number of edges out of the total 
number of edges of M; in this case we can choose the ordering of vertices so that all the relevant indices 
Pi and qi are even, and c = 1. Then we can use the simpler form | det{K^ l )\ = oe det^^, 1 ). 
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The defining property of C(v\,V2) is that it satisfies: KC(yi,V<z) 
function 

1 if V2 = Vi 
otherwise. 

We have the following. 



8 vi (v2)- Here 8 V1 is the delta 



(5ui (V2) 



Lemma 8 T/ie function C is symmetric: C{vi,V2) — C(v2,Vi). We have C(vi,V2) — whenever Vi 
and V2 are both black or both white. If Vi is white, the coupling function C (1)1,1)2) is discrete analytic as 
a function of 1)2, with a pole at V\. 

Proof. Since we already have KC(vx,V2) — 8 V1 (V2), it suffices to show that C(vi,V2) is real when 
V2 — Vi = (1,0) mod 2, pure imaginary when V2 — Vi = (0, 1) mod 2 and zero in the remaining cases. 

If we order the vertices of M in such a way that all the Wo are first, then W\ then Bo and then Si, 
then the matrix K in this basis has the form 



/ 



K 











iK\ 



iK 3 





iK 2 
K 4 







where K\, K2, K3, K4 are real matrices. The conjugate of the above matrix by the matrix 



(I 





V 




il 






\ 




il ) 



is real. Hence the inverse of K has the same form as K. This completes the proof. 



□ 



See Figure M for (part of) an example. 

Since C(vi,V2) = when vi,i>2 are both black or both white, and C{v\,V2) — C(vi,V2), we will 
almost always take the first argument of C to be a white vertex and the second to be black. 



4.1 Boundary conditions for the coupling function 

A discrete analytic function is determined by its boundary values, since its real and imaginary parts are 
harmonic. In this section we describe the behavior of C(vi,V2) for V2 on the boundary of M. 

Assume that vi £ Wo- By Lemma pi G(vi,V2) is real when V2 G Bo(P) and pure imaginary when 
V2 S Bi(P) (and zero when V2 £ Wo U Wi). Let Y be the set of vertices in Bo adjacent to (a white vertex 
of) M but not in M (that is, at distance 1 from a vertex of M). Let Bo'(P) be the graph whose vertices 
are Bo(P) U Y, and whose edges connect every pair of vertices of distance 2, provided that the white 
vertex lying between these two is in M . The set Y is the set of boundary vertices of Bo'(P). Let V 
be the set of exposed vertices di, . . . ,dk (recall that they are all in Bo). See Fig. H for an example of a 
graph Bo'(P). 

Lemma 9 For a fixed vi £ Wo, consider C(vi,V2) as a function of V2. The real part of C(v\,V2), 
extended to be zero on Y and considered as a function on the graph Bq'(P), has the following properties: 

1. it is harmonic at all vertices in Bo(P) \ (V U {vi + 1, Vi — 1}). 

2. AReC{v 1 ,v 1 ±l) = ±1, 

3. its harmonic conjugate is single-valued. 
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Figure 5: Example of the graph Bq'(P) for the polyomino P of Figure @ (P is in dashed lines). The smaller 
gray dots are vertices in Y; the black vertex is the exposed vertex. 

If rather Vi £ W\ then the imaginary part of C(vi, V2), extended to be zero on Y and considered as a 
function on Bo (P ), has the following properties: 

1. it is harmonic at all vertices in Bo(-P) \ (V U {vi + i, Vi — i}). 

2. MmC (vi,vi ±i) = Tl, 

3. its harmonic conjugate is single-valued. 

Proof. The first two properties in both cases follow from 

AC(vi, ■) — K KC(vi, ■) = K*S V1 — S vi+ i — S vi -i — iS vl+i + i5 vi -i- 

This equation is valid at every vertex of Bo(P) except the exposed vertices (which do not have 4 neigh- 
bors). The third property in each case follows by definition, since ImC(«i, •) is the harmonic conjuate of 
ReC(i>i, •) and — ReC(«i, ■) is the harmonic conjugate of ImC(vi, ■). □ 

We will see later that KeC'(vi,V2),lmC(vi,V2) are respectively the unique functions with the above 
properties. As a consequence we will be able to use some general theorems about harmonic functions to 
reach conclusions about the coupling function. 

The conditions in Lemma blare particularly simple because we started with a Temperleyan polyomino. 
For a polyomino with different boundary conditions, the corresponding boundary conditions for the 
coupling function can be quite complicated: see section Kl 



5 Asymptotic values of the coupling function 

Here we will show that, as e tends to 0, the scaled discrete analytic function -C(vi, ■) converges to a 
pair of complex-analytic functions Fo,F\ (Fq when vi £ Wo and Fi when vi € W\) which transform 
analytically (see Proposition ha) under conformal mappings of the domain U. 
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We first study what happens when the polyomino P is the whole plane, since as we will see, for any 
region U the leading term in C(v\,V2) equals Co{vi,V2), the coupling function on the plane (as long as 
vi is not too close to the boundary of U). 



5.1 On the plane 

In [|13| we gave an explicit formula for the coupling function on Z 2 . This was shown to be the limit 
as n —* oo of the coupling function on the 2n x 2n square, centered at the origin. In that paper we 
used different weights for the Kasteleyn matrix: 1 on all horizontal edges and i on all vertical edges. 
The present calculation is straightforward using the same methods (in fact the result is identical after 
changing the sign on alternating vertices of -Bo an d B\) and yields the following. 

„ ^2 



Proposition 10 (|13|) Let Co denote the coupling function for the whole 

i(x8-y<j>) 



Then 



C (Q,x + iy) = 



4tt 2 



2isin(6>) + 2sin(0) 



By translation invariance, Co(vi,V2) — Co(0,V2 — Vi) so this theorem describes the entire coupling 
function. In |l3] it is shown how to evaluate explicitly this integral. Figure M shows the first few values of 
Co(0, x + iy) when x + iy is in the positive quadrant. The values in the other quadrants are obtained by 
the symmetry Co(0, iz) — — iCo(0, z), which arises from the corresponding symmetry of the edge weights. 



13 10 , 



4 5_ 
lit" 4 o 



(0,0) 



... 1 1. 
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Jt " 4 :0 



Mi> 



7U' 4 



4 It 
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4 " 71 



13 10^ 

4 " n 



Figure 6: The function Cq(0,x + iy), the coupling function for Z 2 . 

Recall that the origin in Z 2 is a vertex of type Wo- 
Theorem 11 As \z\ — » oo, the coupling function on 1? is asymptotically equal to —, that is 

{ Re^ + OtAj) 2 G B 

0.(0..)= «™i + o{V) .ea. 
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Proof. There is the following relation between Co and the Green's function for the plane. The real 
part of Co is the unique function on Bo(Z 2 ) satisfying A_ReCo = Si — <5_i and tending to at infinity 
(see Lemma fil and recall that Co is the limit of C on square regions centered at the origin) . 

Now the classical Green's function Go(v, w) on 1? satisfies AGo(0, w) = So(w) and for any fixed v, 
Go(0, w) — Go(v, w) — » as w — » oo (see Lemma [12|). As a consequence we have 

ReC (0,™) = G o (0, H^l) - G (0, H±±), 

where on the right we used coordinates on Bo(Z 2 ) which has index 4 in Z 2 . 
Using Lemma [12^ we have 

ReC (0,™) = Go(0,^i)-G (0,^-i) 

= 2^ V 10 * 1- a" | - l0B| — a" ^ - | -° ( |^F ) 

= ^Relog(^±i) + 0(^) 

2-7T W> — 1 |w| 2 

1 2 1 

= —Re 7 + 0(^) 

lit w — 1 |wp 

= Re— +0(^) 
where we used log(l + z) = z + 0(|z| 2 ). A similar argument holds for the imaginary part. D 

Lemma 12 ([jig]) For the Green's function Go on Bo(Z 2 ) we have 

Go(0,u) = ~logH + co + 0( r ^) (3) 

27T \V\ 

for a constant Co- 

Note that Stohr's Laplacian is —1/4 times ours, so his Green's function is —4 times that in (H). 

5.2 The half-plane 

For later use we will need to compute the coupling function on a half- plane. Let {Pn} be a sequence of 
Temperleyan polyominos in the upper half plane H — {x + iy £ Z 2 | y > 0}, such that P n contains the 
rectangle [— n, n] x [1, n], and the base point do of P n is outside this rectangle. Then (as we will show in 
the proof of Theorem |l4[), for fixed l>i,«2 the coupling function C^ n '(v\,Vi) on P n converges to a limit 
Ch(vi,V2) satisfying the properties below. In particular the uniform measures on the P>"' converge to 
a unique measure fiH- 

Suppose vi 6 Wo- The real part of Ch (vi^Va) satisfies the conditions of Lemma In AReCsr(ui, •) = 
8 vi +i — Sv x -i, ReCjj(wi,a; + iy) — when y = 0, and ReCr/ tends to zero at infinity. There is a unique 
harmonic function with these three properties: the real part of Co(«i, Vi) — Co(W, V2) (note that Vi G Wo 
implies vT £ Wo)- The conjugate harmonic function IraCn is single- valued, and uniquely defined by the 
condition that it tends to zero at infinity; as a consequence we have 

Ch(vi,V2) — Co(vi,V2) — Co(vT, V2) when vi £ Wo- (4) 
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If Vi 6 Wi, on the other hand, it is the imaginary part of Ch(vi,x + iy) which is zero when y = 0. In 
this case there is again a unique harmonic function satisfying the requisite properties: KeCn (i>i,t>2) = 
Re(Co(vi,V2) + Co(VT,V2))- So then 

Ch(vi,V2) — Co{vi,V2) + Co(vT,V2) when vi €Wi. (5) 

There is a big difference between these two cases: from Theorem |llj, in the case vi £ Wo we have 

C H (v 1 ,v 2 ) = -(— ^=)+°(i ' 



TT \V2 — Vl V2 — Vl J |t>2— Vl] 

Vl ~ w +o( T 1 



12- 



7r(v2 — Vl)(l>2 — Vl) \V2 — V\\ 2 

which is O(d), where d is the distance from vi to the boundary. In the case vi £ W\, rather, we have 
C H (vi,v 2 ) = -( — + L_]+0( ' 



TT \ V2 — V i 1)2 — Vl I ' If 2 — Vl I 2 ' 



a*-*-* +o{-^—,) 



TV(V2 — Vl)(v2 — Vl) \V2 — Vi\' 2 ' 

which does not go to zero as Vi approaches the boundary. 

There are similar formulas for the other half-planes with horizontal or vertical boundary. 

5.3 Bounded regions 

One of the main results in this paper is to show that the coupling function on a finite region converges, as 
e tends to zero, to a pair of analytic functions which transform analytically under conformal maps of the 
region. For a fixed region U we can not prove this for all Temperleyan polyominos P e approximating U: 
we require that the approximating P e have a nice behavior in a neighborhood of their exposed vertices. 
This shortcoming is due to our lack of understanding of the asymptotics of the discrete Green's function 
near the boundary of a polyomino. It seems nonetheless reasonable to suspect that this flaw can and will 
be overcome in the near future. 

We will begin at this point to use the metric on eZ 2 rather than Z 2 . That is, we work on polyominos 
in eZ 2 with interior dual graphs having edges of length e. The graphs Bo'(P) have edges of length 2e. 

Let U be a region in C with smooth boundary (or piecewise smooth as previously defined). Let 
Do, • • • , Dk be the boundary components of U, with Do being the outer component. Let d'j be a marked 
point of Dj. Let zi be a point in the interior of U and Z2 be any point of U. 

We define two functions Fo(zi,Z2) and Fi (21,22), whose existence and uniqueness will be shown in 
the proof of Theorem FLU below. For fixed Zi, the function ^0(21, 22) is analytic as a function of 22, has 
a simple pole of residue 1/tt at 22 = z\ and no other poles on U except possibly simple poles at the 
d'j, j > 0. Furthermore it is zero at d' and has real part on the boundary of U. For fixed Zi, the 
function .£1(21,22) is analytic as a function of 22, has a simple pole of residue 1/tt at 22 = 21 and no 
other poles on U except possibly simple poles at the d'j, j > 0. Furthermore it is zero at d' and has 
imaginary part on the boundary of U. 

For each e > sufficiently small, let P e be a Temperleyan polyomino in eZ 2 approximating U in the 
following sense. The boundaries of P t are within 0(e) of the boundaries of U, and except near a corner 
of dll the tangent vector to dU points into the same halfspace as the direction of the corresponding 
edges of dP e . Furthermore assume that the exposed vertices dj of P e are within O(e) of the d'j. Suppose 
further that for a certain 8 — 5(e) > tending to zero sufficiently slowly (see below), in a ^-neighborhood 
of each dj, the boundary of P e is straight (horizontal or vertical). Let M e be the interior dual of P t . Let 
vi be a white vertex and V2 a black vertex of M e . We then have the following result. 
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Theorem 13 Fix any real £ > 0. The coupling function C(vi, V%) on the graph M e satisfies: for vi £ Wo 
and vi, «2 not within £ of the boundary of M t , 

11 

-C(V1,V2) — -Co(vi,V2) + F (V1,W2) +o(l), 

where F * is defined by the condition that Fo(z\,Z2) — , _ z ■, + F * (zi, Z2), with Fo as above, and Co is 
the coupling function on el? . 
If vi £ W\ , rather, then 

-C(vi,V2) - -Co(vi,V2) + Fi(vi,V2) +o(l), 

where F* is defined by the condition that Fi(zi,Z2) — , 1 _ z ; +^(21,22), with F\ as above. 

The equality in the theorem should be interpreted as saying: when vi £ Wo and V2 £ -Bo then 
C (1)1,1)2) equals the real part of the right-hand side; and when vi £ Wo and V2 £ Bi then C(vi,V2) 
equals i times the imaginary part of the right-hand side. Similarly for vi £ W\ and V2 £ -Bo, then 
C(vi,V2) equals i times the imaginary part of the right-hand side; when vi £ W\ and V2 £ Bi then 
C(vi, V2) equals the real part of the right-hand side. 

When vi and V2 are far apart (not within o(l)) then we can replace |Co(t>i,t>2) with , \ v . + o(l) 
and so the statement is simply 

-C(vi,v 2 ) = F :i (vuV2) + o(l) 

where j — or 1 as the case may be. 

Proof. Let Us be equal to U except in a 25- neighborhood of the d'j , and such that Us is flat and 
horizontal or vertical in a S- neighborhood of the d'j . We will first prove the theorem for Us for any fixed 
5>0. 

We will do only the case Vi £ Wo- The case vi £ Wi is identical using the imaginary part of C rather 
than the real part of C below. 

Let G(wi, 11)2) be the Green's function on Bo'(P<!) (recall the construction of Bo'(Pt) from section [4. l[ ), 
that is, the function which satisfies AG(w\, W2) = S W1 (1^2) and G(wi, (02) = when W2 £Yl)V\ {w±}. 

The function R&C(vi,V2), considered as a function of 1)2, is a linear combination of the Green's 
functions G(vi ± e, V2) and G(dj,V2) for j — 1, . . . , k since it is harmonic off of these vertices. In fact 

since 

fc 

AReC(ui, •) = 5 V1+C - 5 Vl - e + ^J ctj5 dj 

3=1 



ReC(vi,V2) =G(vi +e,v 2 )-G(vi - e,v 2 ) +^ajG(dj,V2). (6) 

3=1 

By Corollary llSI below, the rescaled Green's function -G(dj,l>2) (considered as a function of 1)2) 
converges away from dj to a continuous harmonic function with a logarithmic singularity at cL and 
boundary values 0. (This is the place where we need Us rather than U .) Similarly by Lemma llT] the 
difference ^(G(vi + 6,1)2) — G(i)i — €,1)2)) converges. It remains to show that the coefficients a.j in (n) 
converge as e -—> 0. Note that if U is simply connected then k = and we are done. 

For general U, the right hand side of (n) automatically satisfies the conditions (1) and (2) of Lemma 
^defining the coupling function, but the Green's functions G(vi,V2) do not in general have single-valued 
harmonic conjugate. It is necessary to choose the olj so that the harmonic conjugate of the right-hand 
side of (|q) is single-valued. We show that in fact the Oj are uniquely determined by this property. 
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We will use the language of electrical networks, see e.g. |8|. Consider the graph Bo'(P) e to be a 
resistor network with resistances 1 on each edge. The function G(vi,V2) is the potential at V2 when one 
unit of current flows into the network at vi and the boundary FuVis held at potential 0. The aj must 
be chosen so that, when currents aj flow into the network at dj, and current ±1 flows into the network 
at vi ± e, and the boundary is held at potential 0, then the net amount of current exiting each boundary 
component Dj is zero. For, the harmonic conjugate is the integral of the current flow: the integral of the 
current crossing a closed curve surrounding Dj is if and only if the harmonic conjugate is single-valued 
around that curve. 

We claim that given any k + 1 real numbers co,Ci, . . . , c^ such that Co + • • • + Ck = 0, there exists 
a unique choice of reals an,..., a* such that, when currents aj flow into the network at dj, and the 
boundary is held at potential 0, the net current flow out of each boundary component Dj is Cj. This 
will then determine the otj , because letting co , . . . , c* be the current flow out of the boundaries from the 
function G(v\ + e, V2) — G(v\ — e, V2) (we mean, when 1 unit of current flows in at vi + e and 1 flows out 
at v\ — e), we must choose the unique aj to exactly cancel this flow. 

To prove the claim, note that the map $ : R ft — > R fc which gives the outgoing currents ci, . . . , c*. 
(and therefore Co = — ci — • • • — Ch as well) as a function of CKi, . . - , a* is linear (this is the principle of 
superposition). It suffices to show that the determinant of <E> is nonzero. 

However on each column of the matrix of "J> (in the basis {ci, . . . , c^} and {ai, . . . , a^}) the diagonal 
entry is the only negative entry: G(dj,Vo,) induces a positive net current flow out of each boundary 
component except the component Dj which contains dj, since G{dj,V2) is a positive harmonic function. 
Furthermore the diagonal entry in "J> is larger than the absolute value of the sum of the other entries 
in that column, since a nonzero amount of current flows out of Do, i.e. Co > (and the total inflowing 
current equals the total outflowing current). This implies that det$ 7^ (see Lemma |l6| b elow) . 

Now as e tends to 0, the rescaled Green's function -G(dj,V2) converges (Corollary \E)\). This implies 
that the entries of the matrix of $ converge: the pointwise convergence of a sequence of harmonic 
functions implies convergence of their derivatives (even in the discrete case), due to Poisson's formula: 
the derivative at a point is determined by integrating the values of the function on a neighborhood of that 
point against (the derivative of) the Poisson kernel. By integrating the derivative we get convergence of 
the net current flow out of each boundary. Furthermore the amount of current out of Do due to -G(dj, •) 
is bounded from below. This implies that det $ is bounded away from (Lemma nil). Since the difference 
in Green's functions -G(vi + e, ■) — -G(vi — e, ■) also converges (Lemma [17|), the net current out of Dj 
from ~G(t>i + e, •) — ~G(vi — e, •) converges. Therefore the aj converge as well. We conclude that ReC 
converges. 

The C°-convergence of ReC implies convergence of its derivatives and so by integrating we get local 
convergence of ImC as well. By uniqueness of the harmonic conjugate (up to an additive constant) we 
have that ImC converges (the constant is determined by the fact that it is zero at do)- 

In conclusion when vi € Wo, -C(vi,Vz) converges to an analytic function (of V2) with all the properties 
of the function Fq . Furthermore the proof shows that there is a unique function with these properties. 
When Vi € W\ then C(vi,t>2) converges to Fi which is also unique. 

When \v2 — Vi\ = o(l), the main contribution to C{vi,Vi) is from G(vi + 6,1)2) — G[vi — £,^2); the 
unrescaled Green's functions ajG{dj,V2) contribute at most o(l). Since G(i>i + e,V2) — G{v\ — e, V2) — 
Go(vi + e, vs)— Go(vi— e, vz)+o(l) (see the proof of Lemma |l7|) , we conclude that G(vi,V2] = Co(fi, 1*2) + 
o(l). This gives the "local" term in the statement. 

The above holds for Us for any S > 0. It remains to see that when 8 — > the functions F , F± 
on Us converge to Fo, F\ on U. This follows from Proposition ha below, and the fact that the Riemann 
map from Us to U converges (if appropriately normalized) to the identity mapping. Therefore the result 
holds for U as long as 5 — > sufficiently slowly. □ 

A similar result holds when v± is close to a flat boundary of P t . Here is the statement when it is close 
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to a flat horizontal boundary. This is the only case we will need later. 

Theorem 14 Fix S > 0. Let z\ be a point on the boundary of U such that the boundary is flat and 
horizontal in a 5 -neighborhood of z\. Let v\ £ Wo be a point within 0(e) of Zi and V2 a black vertex. 
The coupling function C(vi,V2) satisfies 

1 1 

-C(vi,v 2 ) — -C H (vi,v 2 ) +o(1) 

where Ch is the coupling function defined in M) for appropriate half-plane in el? . If rather v\ £ W\ 
then 

-C(vi,v 2 ) = -Ch(vi,V2)+F**(zi,v 2 ) + o(1), 

where F" is defined by the condition that Fi(zi, 22) = , 2 _ z -, + F**(zi, 22) and F\ is as before. 

Proof. We use the notation of the previous proof. If vi £ Wo, then by (0), the function -(Go(vi + 
e, v 2 ) — Go(vi — e, v 2 )) is already o(l) for v 2 near the boundary of U except at the point 21. Therefore the 
ctj will all tend to as well. The result follows if we define Ch(v\, v 2 ) = Go(vi + e, v 2 ) — Go(vi — e, v 2 ). 

On the other hand if vi € W\, then by (H), the function -(Go(vi + ie,v 2 ) + Go(vi — ie, v 2 )) has two 
poles (each of residue 1/tt) within o(l) of Vi ■ The remainder of the proof is similar to that of the previous 
theorem. □ 

Again note that when v 2 and vi are not close, in case vi € Wo we have -C(vi, V2) — Fo(vi,v 2 ) + o(l) — 
o(l) and when vi 6 W\ we have \C(y\,v 2 ) = F\(v\,vi) +o(l). 

The functions Fo, Fi depend only on the conformal type of the domain U in the following sense. Let 
F+ = Fo + Fi and F- = F - Fi. 

Proposition 15 The function F + (z 1 ,z 2 ) is analytic as a function of both variables. The function 
F-(z\,z 2 ) is analytic as a function of z 2 and anti-analytic as a function of z\. If V is another do- 
main with smooth boundary and if f: U — > V is a bijective complex analytic map sending the marked 
points on U to those ofV, and if FJ{_ ,F_ are the functions defined as above for the region V then 

F?(v,w) = f'(v)FX(f(v),fH) (7) 

FH(v,w) = f'(v)FY(f(v),f(w)). (8) 

Proof. We already know that F+,F- are analytic in the second variable. Going back to the coupling 
function, for a fixed black vertex v 2 not adjacent to v\ we have 

—C(vi + e, v 2 ) + C(vi — e, v 2 ) — iC(vi + ie, v 2 ) + iC(vi — ie, v 2 ) = 0. 

If v 2 6 Bo and v\ + e £ Wo this gives in the limit (using Theorem f!3| ) 

-d xl ReF (vi,v 2 ) + d yi ImFi(vi,v 2 ) - 

and if v 2 £ B\ and Vi + e £ Wo this gives 

-d xl ImFo(vi,v 2 ) - d yi ReFi(vi,v 2 ) = 0. 

These can be combined into a single complex equation 

-d xl Fo(vi,v 2 ) - id yi Fi(vi,v 2 ) = 0. 

Similarly if vi + e £ Bi this gives 

—d xl Fi(yi,V2) - id yi F (vi,v 2 ) = 0. 
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Summing these gives dzr(Fo + Fx) — and taking their difference and conjugating gives <9 Z1 (Fo —Fi) = 0. 
This proves the first two statements. 

As a function of Zi, the function Fq (/(zi), / (#2)) has all the properties of Fq except that the residue 
at Z2 = 21 is —777- i ■ Similarly the function F^(f(zx), f{z2)) has all the properties of Fl 7 except that 
the residue at 22 = Zx is — yrV- ; ■ So letting a.,/3 be the real and imaginary parts of /'(zi) we have that 

a(z 1 )F^(f(z 1 ),f(z 2 )) + if3(z 1 )F^(f(z 1 ),f(z 2 )) 

has residue ^»T Z I = ^ at 22 = 2i, and all the other properties of Fq , and so must equal Fq since 
i^o 7 is unique. A similar argument shows that 

i(5(z 1 )F^(f(z 1 ),f(z2)) + a(z 1 )F^lf(z 1 )J(z 2 )) = Fl J . 

The equations for F+ and F- follow. □ 

As an example, on the upper half plane we have from m) and (0) that 

Fo{zi,Z 2 ) = —. r -. =r, 

iv(z2 — Zi) ir{Z2 — zx) 

and 

Fi(zi,z a ) = -7 r + - 



7r(z 2 — Zx) Tv(z'2 — Zl)' 

These functions vanish at oo, which can be thought of as the location of do. In particular Fj r (zx,Z2) — 

—7- 2 _ z •. , which is analytic in both variables, and F-(zi,Za) — r- 2 _ E ) , which is analytic in 22 and 

antianalytic in zx- 

Let U be the upper half plane with do located at (that is, a square of type B\ is removed near 
the origin). We can compute Fq ,F± for this new region U by using the above transformation rules. A 
conformal isomorphism from the upper half plane to itself which takes to 00 is f(z) = — 1/2. 

Since f'(zx) = z^ 2 we have 

F?(zx,z 2 ) = 1 2 



4 *(/(■») - /(*0) 

2z 2 



TYZl(z2 — Zl) 

Any other choice of f(z) would give the same result. The function FY. is obtained similarly. 
Lemma 16 Suppose 5 > 0. If Q is an n x n matrix Q — (qij) and for all i, 



3, 3^i 



then detQ > 5" > 0. 



Proof. Gaussian elimination using rows preserves this property: if for each j we multiply the first row 
by Qjx/qxx and subtract it from the jth row, the first column of the new matrix is all except for the 
first entry qn, and the remaining n — 1 x n — 1 submatrix still has the property in the statement. For 
example the first column of the submatrix is 

<?12 §12 912 
1722 921, <732 ?31, • ■ ■ , <?n2 QnX 

qxx qxx qxx 
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and 



912 j- . I I <J12 

1722 g2i — o > \qi2\ + \qsa\ + ■ ■ ■ + \q n 2\ g2i 

9n qn 

^ f, I I<?i2 1 ■ |gai 



= |?32 H + 9n2| + |gi2| 



gn 

gn - [ggij 
?11 



> |t?32 1 H 1- |qn.2 \ + ^-(S + \q 3 i +--- + \qm\) 

qn 

■^ I <?12 1 , ,, 012 I 

> 932 931 1 H h |<7n2 ?nl - 

□ 

Recall that the continuous Green's function on a region U is the real-valued function gu satisfying 
Agu{zi, Z2) — 8 Z1 (22), and which is zero when 22 is on the domain boundary (here S Z1 is the continuous 
delta-function, and A = —■£?-* — ■£-*). 

or oy ' 

Lemma 17 Let Z\ = x\ + iyi be a point in the interior ofU, and let z z G U, 22 7^ Z\. Let vi be a vertex 
o/Bo'(P e ) within O(e) of z\, and let v% be a vertex o/Bo'(P e ) within 0(e) of z%. Then the difference of 
(rescaled) Green's functions -G(vi + e, V2) — -G(v% — e,V2) converges to 2d Xl gu{zi,Z2) ■ 

Proof. Let H(vi,V2) = -(G(vi + e, V%) — G(vi — e,V2))- From Theorem [ill on the plane eZ 2 we 
have 

Ho(vi,V2) = -(G (vi + e,v 2 ) - G {vi - e,v 2 )) = Re— r- + 0( -. 7-). 

e ir(V2 — Vi) \V2 — Vip 

The function H(vi,V2) — Ho(vi,V2) is harmonic (as a function of V2) on all of Bo' (P e ) (including Vi ± e) 
and has bounded boundary values, since H (vi,v 2 ) is O(l) on the boundary of B '(P E ) and H(vi, v 2 ) is 
zero there. Let g be the continuous harmonic function which has boundary values equal to the boundary 
values of the limit 

lim H(vi, v 2 ) - H (v 1 ,v 2 ). 

Since these boundary values are continuous in the limit, g exists and is unique. Note that the boundary 
values of H — Ho are within O(e) of the limiting values (Theorem |ll|). 

Restrict jtoa function on the vertices of Bo'(P e ). The discrete Laplacian of g at a vertex v £ Bo'(Ps) 
is: 

A e g(vi,v) =4g(v) -g(v + e) - g(v - e) - g(v - ie) -g(v + ie) 

and when e is small we can approximate this using the Taylor expansion of the smooth function g, 
yielding 



• ' fd 4 q(v) <9 4 <7(u)\ , , 



e^_ 

24 V dx 4 ' dy 4 

Therefore H (1)1,1)2) — Ho(vi,V2) — g(vi,V2) has discrete Laplacian which is 0(e 4 ) on Bo'(P e ), and 
the boundary values are 0(t). A standard argument now shows that H — Ho is close to g: the function 
x + iy 1— > x 2 has discrete Laplacian which is a constant; choose constants B2,B% sufficiently large so that 

A E (B 2 e 4 (Re(v2)) 2 + H(vi,v 2 ) - Ho(vi,v 2 ) - g(vi,v 2 )) >0 

and 

A, (P 3 e 4 (Re(« 2 )) 2 -H(vi, v 2 ) + H (vi,v 2 ) + g(v u v 2 )) > 
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on Bo'(-Pe). By the maximum principle for superharmonic functions, these functions must take their 
maximum value on the boundary of the domain Bo'(Pe). Since H(vi,V2) — Ho(vi,V2) — g(vi,V2) = 0(e) 
on the boundary of Bo'(P E ), we conclude that 

\H(vi,v 2 ) - H (vi,v 2 ) - g(v 1 ,v 2 )\ = O(e). 

Therefore H(vi,V2) converges to the function Re — j- 1 _ v . + g(vi, V2) which has boundary values and a 
single "pole" of residue 1/tt at vi. This is 2 times the a; 1 -derivative of the continuous Green's function. 

□ 

A similar result holds for the j/i-derivative of gu, yielding: 

Corollary 18 Recall the definitions of the functions Fo,Fi,F+,F- from Theorem |p] and Proposition 
\lq. Letting Z\ — Xi + iyi, we have 

2dgu(zi,z 2 ) — Fo(zi,Z2)dxi + Fi(zi,Z2)dyi = -F + {z\, z 2 )dz\ + -F- (zi,za)dzi 

where the exterior differentiation dgu is with respect to the first variable. 

When z\ € dU the proof of Lemma hj implies the convergence of the Green's function as well. 

Corollary 19 Let 5 > 0. If Zi is on the boundary of U , and the boundary of both U and P E ts straight 
and horizontal in a 8 -neighborhood of z\, then for v\ within O(e) of zi, 

-G(V!,V2) = 3(7(21,22) +0(1). 

Proof. Reflect Bn'(P e ) across the boundary edge near z\ (the edge consisting of vertices in V) to get 
a graph Bo"(P e ). Glue Bo'(P £ ) and Bo"(P e ) along their common edge in a 5-neighborhood of Vi. A 
harmonic function / on Bo'(P e ) which is zero on the boundary extends to a harmonic function on this 
glued graph by setting /(«') = —f(v) when v' is the reflection of v. In other words the Green's function 
G(vi, V2) on Bo' (P £ ) is the difference of two Green's functions on Bo'(Pe) U Bo"(P e ); one centered at vi 
and one centered at v[ . 

On the glued graph Bo'(P e ) UBo"(Pe), the vertices vi,v[ are at distance at least S from the boundary 
<9(Bo'(P c ) U Bo"(P e )), but only distance O(e) from each other. The argument of Lemma |l7] can then be 
applied in this case, replacing H(vi,V2) by -(G(vi,V2) — G(vT, ^2))- □ 

A similar result holds when the boundary is vertical. 

6 Conformal invariance of heights 
6.1 Proof of Theorem [j] 

Let U be a region in C with boundary which is piecewise smooth as previously defined. Let dj be a point 
on the j'-th boundary component Dj of U. Let e'j 7^ dj be another point of Dj, which is not at a corner 
of the boundary. 



Let P e be a Temperleyan polyomino approximating U in the sense of section 5. J , with the additional 
constraint of having horizontal boundary in a neighborhood of each e'j , and so that the interior of U is 
locally below each ej . We show that the distribution of the heights of the boundary components of P e is 
conformally invariant. 

Let ej be a vertex on the boundary of P e near e'j . We assume for simplicity that each ej has the same 
parity (its coordinates have the same parity) as en. For definiteness we suppose the lattice square whose 
lower left corner is Cj is of type B\ for each j. 
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Let hj be the random variable giving the height of e-j for a random tiling of P e assuming the height 
of eo is zero. Let hj be the mean value of hj . 

We will show that for integers ni,n,2, ■ ■ ■ ,Uk > 0, the moment 

E((fti - M ni (Aa - M" 2 • • • (h k - h k ) n ") (9) 

is conformally invariant. Let K — ill H — • + rife. The precise value of the moment (]9t) is as follows. 

Proposition 20 Let {ji}i£[i,K] be a collection of pairwise disjoint paths in U , such that for each j G [1, k] 
there are ru paths runnning from the outer boundary to the jth boundary component. Then as e — > the 
moment ([^) converges to 



E 



ei ■ 



■£k 



.,e K e{±l} 



det (f £ 

i,jE[l,K]\ ■ 



,(Zi 



dz\ 



(si) 



■ dz 



(e fc ) 



(10) 



where dz- = dzj and dz- = dz7, and 
j j j j 



F e 



if i = j 

F+{ Zl ,z 3 ) if ( £l , £i ) = (l,l) 

F-j^Zj) if (£„£,) = (-1,1) 

f-(^,^) if (ei, £,-) = (1,-1) 

L F+( Zl , 2j ) if (Ei,E i ) = (-1,-1). 

Note that in each of the 2 K multiple integrals in (|l0|), the integrand I is conformally invariant, in the 
sense that 



I(z)dz ■■ 



/(/(z))dz. 



/(7) 



This follows because of the transformation rules (vti and the fact that each integrand is analytic or 
antianalytic in Zi according to e; = ±1. Ther efor e the moment (H) is conformally invariant. 

An example calculation is done in section p.3| . 

By H, section 30], there is a unique probability distribution with these moments on condition that 
the moment generating function 



H(ti 






m(m, . . . ^fe)*" 1 



mi 



■n k \ 



has nonzero radius of convergence around the origin (here m(ni, . . . ,n k ) is a shorthand for (RSI)). This 
convergence is shown in Lemma E3, below. We can then conclude that the probability distribution with 
these moments is conformally invariant, and by M, Theorem 30.2] that this distribution is the limit of 
the distributions for finite e. This will complete the proof of Theorem 111. 
Proof of Proposition BQ. 



For each e sufficiently small and for each j £ [1, k] let 'y \ 



(0 



„( £ ) 



be pairwise disjoint lattice paths 



(which are also disjoint for distinct js) in P t which start on the flat boundary near eo and end on the 
flat boundary near ej . We require that each straight edge of 7^ have even length (by this we mean, a 
length which is an even multiple of e) . This is possible by our choice of parities for eo and ej . 

In a given tiling the height change on 7^ equals 4(A, S — Bj S ), where Aj B is the number of dominos 
crossing 7^ with the black square on the right and Bj S is the number of dominos crossing *ffj with the 
black square on the left. To see this, note that if 7;*' does not cross any dominos, the height change is 
0: the straight edges have even length so the height change along them is zero. Then, for each domino 
crossed by JjJ, the height difference changes along that edge from — 1 to +3 if the domino has black 
square on the right, and from +1 to —3 if the black square is on the left. 
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Since hj = 4(A,- a — Bj s ) for each s, the moment (H) is equal to 

4 K E((A n - B n - A 1± + Bn) ■ ■ ■ (A knk - B k „ k - A k „ k + B fen j) (11) 

where K — m + • • • + rife. 

The remainder of the proof involves expanding this out, cancelling various terms and then recombining 
in the right way. 

For notational simplicity we renumber the paths 7^ from 1 to K. Similarly change indices of Aj B , Bj B 



to values in [1,-K]. For j £ [1,-K] let ctjt be the t-th possible domino of 7 crossing 7 whose black 



» 



» 



(<0 



square is right of 7 . Similarly let (3jt be the t-th possible domino crossing 7 whose black square is 
on the left. Let oijt,f3jt also denote the indicator functions of the presence of these edges/dominos. Then 

^-■B*=£>t-x;&*. (12) 



Let (u!j s ,bjs) be the white and black squares, respectively, of the domino Oij s and (w'j s , bj s ) be the 
white and black squares of the domino /3j s . 

Since the straight edges in the path 7] have even length, we can pair the ayt dominos with adjacent 
j3j t i dominos which are parallel to a, t . It is then convenient to write 



Bi - Aj + B 3 =Y. (^' - & i* - &* + M 



where a.jt and Pjt are paired. Equation (U/w is now 



4 2_^ ^[(Cxiti ~ Otlti — Pit! + 01tt) ■■■ (aiKt K — OLKt K — 0Kt K + 0Kt K )) , 
t± ; ...,tg 

where the sums are over all pairs {(Xitx-tPux) °f 7i 1 («2t 2 i/^2t 2 ) °f I2 an< i so on - 
Lemma 21 Let e% = (u)i, bi) for i — 1, . . . , n be a set of n disjoint edges; then 

( C(wi,b 2 ) ... C(wi,b n ) \ 



E((ei - ei) • • • (e n - e n )) = a E det 



C(w2,bi 







C(w„-l,b n ) 

C(w n ,b n -i) 



V C(w»,bi) 

where (using the convention after Theorem it) oe is the product of the edge weights of the e 
Proof. This follows from Theorem M, induction on n and the fact that 



ai2 

CJ21 0122 
a n 1 



<lln 



an 012 

0,21 0.-22 



Bill 



a 1 r, 



an 
a22 

a n2 





(I 2n 



(13) 



a 



Now expand the summand of ( |l3| ) into 2 A terms 

E((a ltl - a ltl ) . . . (a KtK - a K t K )) + ... + {-l) K E({f3 ltl - Ul ) . . . {(3 K t K - P~Kt K )). (14) 
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By Lemma Ell, each term is a certain quantity cie times the determinant of a K x K matrix whose entries 
are given by the coupling function connecting black squares of the dominos a s t B , f3st s with white squares 
of the other dominos. Since each '(3' edge has weight of the opposite s ign as the 'a' edge to which it is 
paired, the signs in Mm cancel with the sign changes in the cle and so (TL4J) is equal to the sum of all 2 K 
determinants, times the product cle of the edge weights of the first determinant. 
Consider the first term in ( |14[ ) 

E((a ltl - Qi tl ) . . . [pL K t K - otKt K ))- (15) 

Recall that (wj s , bj s ) = etj a and (w'j s , bj s )_= f3j 3 - Fix a choice of indices s = 
can drop the second subscripts. By Lemma glj, equation (|15| ) is then equal to 



C(wa,6i) 

C(wi,b 2 ) 



C(w K ,bi) 



C(w K ,b K -i) 
C(w K -i,b K ) 



(16) 



C(wi,b K ) 
A typical term in the expansion of (K&) is 

assgn(cr)C , (tui,& .(i))C(tu 2 ,& £ r(2))- ■ -Ciwicb^K)) 

where o has no fixed points. 

Let us first assume that a is a A'-cycle; reorder the indices so that ( |17| ) becomes 

a E sgn(a)C(wi,b 2 )C(w 2 ,b 3 ) ■ ■ -C(w K ,bi). 



(17) 



(18) 



To expand this out, define variables r; = ±1 according to whether Wi £ Wo or Wi € W\, and s; = ±1 
according to whether bi £ Bo or bi G B\ . If we assume that neither wi or 62 is close to the boundary, we 
can then write (see Theorem n3l and the remarks immediately after its statement) 



— — —F {w 1 ,b 2 ) H 



r-i 



Fi(wi,b 2 ) +o(e) 



6(^1,62) = e I 7. «ImH Re 

= -(F + (wi,b 2 ) +r 1 F-(w 1 ,b 2 ) + s 2 F-(w 1 ,b 2 ) + riS2F+(wi, b 2 )) + o(e). 

For each fixed ^ > 0, when neither of wi , b 2 are within £ of the boundary, this approximation holds 
for sufficiently small e. When one or both of wi,b 2 are within £ of the boundary, we only need to know 
that -C(u)i, b 2 ) is bounded by some constant independent of e and £. Then in the sum ( |13[ ) (and in the 
integral (hoi) we can ignore all terms in which some u>; or bj is within £ of the boundary, as these will 
contribute at most 0(£). The boundedness of \C{w\,b 2 ) follows from the convergence of the discrete 
Green's function as in Theorem |l3| 

We can now write (Il8|) as 



e a E sgn((j)\(F + (w 1 ,b 2 ) +r 1 F^(w 1 ,b 2 ) + s 2 F~(w 1 ,b 2 ) + ns 2 F+(w 1 , b 2 )) ■ ■ ■ 

(F + (w K ,b-i) +r K F-(w K ,bi) + siF-(w K ,bi) + r K siF+(w K , bi)) ) + o(e A ) 



(19) 



We obtain a similar expression if we replace (vii,bi) by (w[,b'i), except that the signs of r\ and s\ are 
reversed. In particular if we sum up over all 2 v choices of «_, and (3j (as we need to do to obtain (ll 
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we get 2 K times the sum of those terms in (|19J) which have rt to the same power (1 or 0) as Si, for each 
i. This sum can therefore be written as an error o(e K ) plus 

2~ K e K sgn(a)a E ^ {ris 1 ) < - 1 ~ £l)/2 ■ ■ ■ (r K s x ) (1 " Eir)/2 F £1 , E2 {z%, z 2 )F E2>E3 (z 2 , Z 3 ) ■ ■ ■ F EKtE1 (z K , Zi), 

ei,...,e K e{-l,l} 

(20) 

where F ei>e (z%, Zj) is as defined in Proposition [20[ 

Now in view of replacing the sum ( |l4|) by an integral when e is small, we can replace e by a certain 
phase times \dzj or \dz~j. When the path 7-,- is going east (horizontal and to the right), we have 
2e = dxj = dzj = dzj, and the edge of type a has weight —i, because its upper vertex is white and lower 
vertex black (recall that edges of type a have black vertices on their right). Furthermore tjSj = -1 on 
an east-going path. When the path jj is going west, 2e = — dxj = — dzj — —dzj, the edge of type a 
has weight i, and TjSj — — 1. When the path jj is going north, 2e = dyj = —idzj = idzj, the edge a 
has weight 1, and rjSj = 1. When the path jj is going south, 2e = —dyj — idzj 
has weight — 1, and rjSj — 1. Notice that in each case 2e times the edge weight, times (rjSj) 
—Ejidz^ 3 (recall the definition of dzf from Proposition EOJ). Recalling that cle is the product of the 
edge weights (of the a- type edges), for any choices of the Ej we have 

a E (2e) K (r lSl )^-^ 2 ■ ■ ■ (r K s K ) (1 -^ 2 = H)* £l ■ ■ ■ e K dz^ ■ ■ ■ dz { ^\ 

The sum (EOh is therefore 

4" K (-i) K sgn(a) ^ ei • • ■ s K F eitE2 (z 1 , z 2 )F e2tE3 (z 2 , z 3 ) ■ ■ ■ F eK , ei {zK,z±)dz[ El) ■ ■ ■ dz^ K ■ 

ei,...,eK-6{-l,l} 

(21) 

When a is a product of disjoint cycles we can treat each cycle separately and the result is the 
product of terms like (Ell) involving disjoint sets of indices. Thus when we sum over all (fixed-point 
free) permutations we obtain the formula of the proposition, but without the integral. The factor of 
4~ K cancels with the factor of 4 K in (tl3f), and summing over all pairs gives the integral in (|10[). This 
completes the proof. 



Lemma 22 The moment generating function for the moments \1(\ ) has positive radius of convergence. 

Proof. Letting K = n\ + ■ • ■ + Uk denote the "size" of the moment, it suffices to show that a 
moment of size K is smaller than (cK) K for a constant c. Let 71, ... , jk be the paths of integration in 
(ll0|). We can choose the 7; so that no two are closer than c\/K for some constant ci; indeed, we can 
choose the paths so that the distance between jj and 7,- is at least ci\i — j\/K. Since Fq(zi,z 2 ) and 
Fi(zi, z 2 ) are 0( , \ z , ), in the determinant in ( |l0| ) the ij-entry is at most c 2 K/\i — j\ in absolute value. 
The determinant of a matrix is bounded by the product of the ^2-norms of its rows, and each row of the 
determinant in (hoh has ^2-norm bounded by K(2 + ■§! + ■■■+ fg/21 5 ) = C3 ^ ^ or an °ther constant C3. 
Therefore the sum of the integrals in ( p.0| ) is bounded by cf" K K for a constant Ca- This completes the 
proof. □ 

6.2 The average height. 

Let U C C be a region with piecewise smooth boundary as previously defined. Let h be a point on the 



outer boundary of U, b 7^ do- For each e ^C S let P c approximate U as in section 5.3, but with the 
additional constraint of having horizontal boundary in a ^-neighborhood of b. (We also assume that the 
interior of P e is locally above the boundary at b.) Let z be a point in the interior of U. Let z 6 P c 
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be within O(e) of z and let b' £ dP € be within O(e) of b. We assume that b' and z' are the lower left 
corners of lattice squares of type Bi. Let ^ e ' be a lattice path from b' to z' such that all edges of 7^' 
have even length, and which starts straight and northgoing for a distance at least c5 for some constant 
c. In the notation of the previous section, we have M(h(z)) = 4^E(a s ) — E(/3 S ) where a s ,/3 a are pairs 
of potential dominos crossing the path 7'^ . 

Near the boundary, j^ c ' is northgoing. When a s ,/3 s are within o(l) of the boundary we have 

E(a 8 ) = C(u>s, b 3 ) = Ch{w s , b s ) + 0(e) 

and 

E(/3 S ) = -C(v' s ,w' s ) = -C H («4& s ) + 0(e) 
(note that a B has weight 1 and j3 s has weight —1 when 7' E ' is northgoing). Therefore using Theorem u4 

E(a s -ft) = C H (w a ,b s )+C H (w' s ,b' s ) + 0(e) 

= l+C (W7,b s ) + ^+Co(W s ,b f s ) + 0(e) 

= Co(0,6 s -I(J7) + C*o(0,6' s -<) + 0(e). 

Near the boundary, b s — ®J takes successively values 1 + 2i, 1 + 6i, . . . , 1 + (2 + 4fc)i . . . and b' s — w' s takes 
successively values l + 4i, l + 8i, . . . , l + 4ki .... When we sum over all pairs (w s ,b s ), (w' s ,b' a ) on the path 
7' E ' which are within o(l) of the boundary, the contribution is o(e) plus 

C o (0, 1 + 2i) + C (0, 1 + 4i) + ■ ■ ■ + C o (0, 1 + 2A») + ... = -. 

(This formula can be proved analytically from Proposition (hoj) or more simply by symmetry, noting that 
the average height on the upper half-plane is | given that the height on the boundary alternates between 
and 1.) 

For the terms not near the boundary we have, by Theorem [13, when ^' is northgoing, 

G(Ws, bs) + C{w' B , b' B ) = -+eReFi(zs,z s ) + — + eRcFo(zs,Zs)+o(e) 
= Re(Fl(z a ,z s )e)+o(e) 
where z s is the coordinate of w s and F+ = Fq + F* . Similarly for the other directions of 7 we have 

{Re(— F^(z s , z s )e) + o(e) when 7 is southgoing 
lm(F+(z s , z s )e) + o(e) when 7 is eastgoing 
lm(—F+(z s , Z s )e) + o(e) when 7 is westgoing. 

We can replace e by ^dz s , — \dz ai —^dzs, \dz s respectively according to whether 7 is east-, north-, 
west-, or southgoing. Then all four cases become 

C(w s ,b s ) +C(w' s ,b' s ) = -Im(F+(z s ,z a )dz s ) +o(e). 

The average height is then given by the imaginary part of the integral of 2F+(z, z)dz from 6 to z 
(recall the factor of 4 from the first paragraph of this section), plus |, the constant coming from the 
boundary. This expression does not depend on 8. 

For another region V conformally equivalent to U we have the following. Let / : V — > U be a conformal 
isomorphism. Then F+{zy,Zz) = f'(z\)F+(f(z\), f{zi)) from Proposition |lq. Therefore 



(F^Y(z 1 ,z 2 ) = F|( 2l ,^ 2 )- 



7TIZ2 - «ll 



+ /'(*l) ((*?)*(/(*),/(*,))+- 



n(z 2 - Zl ) J — V V +y -™— *(/(«>) -/(*0) 
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and in the limit as z 2 — * Z\ this is (simplifying using the Taylor expansion of /) 

717' Oi) 



So the average height of z£ V equals the average height of f(z) in (7, plus a term 

(log f'(z))'dz. 



fit) 



This term is — — times the change in total turning (in radians) of the path f(<y) from the path 7. 

This implies that if the path 7 starts at the outer boundary of U , at a point where the tangent 
vector (chosen in the counterclockwise direction) has angle 9 with respect to the horizontal axis (where 
9 £ [0, 27r)), then the average height of a point z £ U is 



1 29 or 

- + — + 21m 

2 7T 



/ F+(zi, zi)dzi. 

■J -y 



Therefore we have 



Theorem 23 Up to an additive constant, the average height of a point z not within o(l) of the boundary 
of U is given by the harmonic function whose boundary values are v°< , where 9(x) is the total turning 
(in radians) of the tangent vector to the boundary on the boundary path going counterclockwise from, d 
to x. 

Note that the boundary values are discontinuous at the point d' Q . 

For example, as noted earlier on the upper half plane when do — 00 the average height of every point 



is ^. When do = 0, rather, then recall that F+(zx,Z2) = — 
height at a point z is (integrating from x = Re(z)) 

E(h(z)) = 



21m 



(z 2 -zi) 



2dzi 

"KZ\ 



So F+(zi,zi) — £j-. The average 



1 4 

- + -Imlog(z/:r) 

Z TV 

- + -arg(2). 

Z TV 



This is the harmonic function with boundary values (on the axis) | to the right of the origin and | = h + 4 
to the left of the origin. Note that on the boundary of the polyomino P e , the height alternates between 
and 1 to the right of the origin and between 4 and 5 to the left of the origin. 



6.3 Example: a second moment computation. 

For a random tiling of the upper half plane with do = 00 we compute the moment K((h(p) — h(p))(h(q) — 
h{q)) for two points p, q. Since h(p) — h(q) = |, this will also give E(h(p)h(q)). 

Let r, s be the vertical projections of p, q, respectively, to the x-axis. Let 71 and 72 be disjoint paths 
running straight from the boundary to p, q, respectively. From Theorem hoi we have 

E((h(p) - h(jp)){h{q) - h(q)) = 



F+ (zi,z 2 ) 
F+(z 2 ,zi) 


dzidz2— / 

•'71 ,72 

dzid22+ / 

■'71,72 


F-(z u z 2 ) 


dz\dz 2 


F-(z2,zi) 


F-{zi,z 2 ) 
F_(z 2 ,2i) 


F+(zi,z 2 ) 


dzTdz2 


F+(z 2 ,zi) 



Conformal invariance of domino tiling 27 

For the upper half- plane we have F+(zi, 22) = , 2 _ z > and F_(zi, 22) = ^i z 2 _— ■, ■ Plugging these in 
gives 

4 f f 1 , , iff 1 _, 

-dz\dZ2 H 5- / / 7 =3 - 5 -d2rd22 + 



(22 -2i) 2 7T 2 J 71 J 72 (22 -2j) 

H « I / 7= -rd2id22 / / 7== — = r—dz~[dz2. 

I" 2 J 71 J 72 (22 - ^l) 2 7T 2 7 / (22 - 2l) 2 



The first of these integrals gives 



4 l 0g ( p ~ g ^ r " s ) 
7r 2 (p — s)(r — q) 



Therefore 



m(p)-Hp))(Hq)-h(q))) = -^f-2Relog [ P ^ S | +2Relog ^ q)( f a \ 

n 2 \ (p-s)(r-q) (p-S)(r-q) 

= -j Re log ( 

7r z \p — q 



7 Trees and winding number. 



A directed spanning tree on a (undirected) graph G is a connected contractible (acyclic) collection 
of edges of G, where each edge has a chosen direction such that each vertex but one has exactly one 
outgoing edge. The single vertex with no outgoing edge is called the root of the tree. If G is a graph with 
boundary, (that is, there is a subset of vertices called the boundary of G), then a directed essential 
spanning forest is a collection of edges of G, each component of which is contractible, where each 
edge has a chosen direction, such that each non-boundary vertex has exactly one outgoing edge, and no 
boundary vertex has an outgoing edge. 

"Temperley's trick" (see H) is a mapping between domino tilings of certain polyominos and directed 
essential spanning forests of associated graphs. In the case P i s a Temperleyan polyomino, the directed 



essential spanning forest is on the graph Bo'(-P) of section 4.1 and the boundary consists of the set Y. 
The forest is defined from a tiling as follows. Each square v in Bo PI P is covered by a domino. The 
white square of this domino lies over an edge of Bo'(-P). This edge is chosen to be the outgoing edge of 
V on the tree on Bo'(-P). See Figure M for the directed essential spanning forest associated to the domino 
tiling of Figure H. 

To see that the essential spanning forest constructed from a tiling has no cycles, it suffices to construct 
the planar dual forest, which is constructed in a similar way from the graph Bi(P) U {do}. In the case 
P is a Temperleyan polyomino, the dual forest is a tree rooted at do (since do is the only possible root). 
Since the dual tree is connected the primal tree has no cycles. 

Conversely, any essential spanning forest on Bo'(P) gives a domino tiling of P, so these systems are 
in bijection. 

The height function of a domino tiling has a nice interpretation for the directed paths in the associated 
spanning tree. To a vertex v in Bo'(P) associate a height which is the average of the heights of the four 
vertices of P adjacent to v. If the outgoing edge of the tree at v points to an adjacent vertex v' , and 
the outgoing edge at v' points to a vertex v", then the height at v' equals the height at v if the three 
vertices v,v',v" are aligned; if the path turns left at v' then the height at v' is one less than the height 
at v; if the path turns right at v then the height at v' is one more than the height at v. 

Therefore the height function along the directed path measures the net turning of the path. 

Proposition 24 Let P be a Temperleyan polyomino with a tiling and let T be the associated essential 
spanning forest. The height change along a directed path 7 in T equals the net turning of the path, that 
is, the number of right turns minus the number of left turns. 
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Figure 7: The directed essential spanning forest (solid arrows) associated to the tiling of Figure 0. The dual 
tree is shown in dotted arrows. 

In particular if 7 is a directed path in T running between dj £ Dj and the outer boundary, the height 
difference between Dj and Do is exactly measured by the winding number of the path 7 (around Dj). 

In Figure g we show the spanning tree associated to a tiling of a Temperleyan annulus in which the 
height difference between the boundaries is 4. The directed path from a vertex adjacent to di to do is 
highlighted. 

8 Other boundary conditions 

There are a number of intuitive ideas in the proof of Theorem hi which are worthwhile exploring. Foremost 
is the interesting link between the height function along a boundary component and the singularities of 
the coupling function. When we introduced the exposed vertices in our polyominos (in order to make 
it tilable) we 'created' poles in the coupling function at those points. There are a number of other, 
equally simple, boundary conditions which give different boundary behavior for the coupling function. 
The most natural seems to be to have all boundary edges have even length. This is natural from the 
point of view of tilings since it is trivial to show that such a region has a tiling. Furthermore the height 
function along such a boundary is particularly simple in this case. However the boundary conditions for 
the coupling function are more difficult: on some boundary edges the real part will be zero and on others 
the imaginary part will be zero. The coupling function will have poles at certain corners and zeros at 
the remaining corners. It seems more difficult to prove the convergence of the coupling function when 
e — > in this case. 

Another potential improvement in the proof would be a more general result (more general than 
Corollary [191) concerning the convergence of the discrete Green's function centered near the boundary of 
a domain. Surprisingly, this problem does not seem to have been considered in the literature. 

Another direction to be explored is the case of regions without boundary. In |l3] we computed a 
formula for the coupling function on a torus. By a recent result of Tesler |lq| higher-genus surfaces can 
be handled by similar methods. 
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